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FOREWORD 


This  is  the  second  semi-annual  technical  report  covering 
research  performed  under  Office  of  Naval  Research  Contract  No. 
N00014-66-C0030,  ARPA  Order  749,  entitled  "Analysis  of  Ionospheric 
Data".  The  contract  period  being  covered  is  May  15,  1966  to 
November  15,  1966.  The  contract  is  under  the  technical  cognizance 
of  Mr.  R.  Gracen  Joiner  of  ONR  and  I.t.  Col.  J.  Hill  of  ARPA. 

During  this  report  period  the  individual  analytical  tasks 
described  in  the  previous  technical  report  which  contribute  to  the 
analysis  of  ionogram  data,  high  frequency  phase  and  amplitude 
data,  and  cosmic  noise  absorption  data  were  continued.  The  use 
of  analyzed  data  in  constructing  an  ionospheric  electron  density 
profile  model  for  selected  data  intervals  was  initiated.  This 
report  mainly  consists  of  a  description  of  the  verification  of 
the  ionospheric  model  for  ionospheric  profile  change  at  dawn  for 
the  days  15  to  19  October  1962.  One  day,  18  October,  is  selected 
to  show  the  detailed  agreement  of  the  ionospheric  electron  density 
profile  obtained  with  true  height  analysis  of  ionogram  data. 

The  author  wishes  to  thank  Mr.  J.  W.  Wrigh'  of  ITSA  who 
furnished  the  true  height  analyses  of  the  ionogram  records  utilized 
in  this  report. 

Respectfully  submitted, 

IIT  RESEARCH  INSTITUTE 


APPROVED  BY: 


riilTam  0.  Davies 
Assistant  Director 
Geophysics  Division 


Robert  D.  Sears, 

Senior  Scientist-Group  Leader 
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ABSTRACT 


A  technique  for  determining  the  electron  density  profile 
in  the  E-region  from  group  path,  phase  path,  amplitude  and  true 
height  analysis  data  is  described.  A  selected  data  period, 

1600  to  2100  GMT  on  15  to  19  October,  is  summarized  in  the  detailed 
reduction  and  analysis  of  the  individual  data  obtained  from  the 
different  measurement  techniques.  One  period  of  detailed 
electron  density  profile  computation  is  presented.  The  period 
1700  to  1800  GMT  on  October  1962  also  included  a  moderately 
intense  solar  flare  so  that  the  excess  ionization  caused  by  the 
flare  in  the  lower  ionosphere  may  be  observed.  The  true  height 
analyses  which  are  developed  are  extensions  of  the  true  height 
analyses  of  the  F-region  and  serve  to  extend  conventional  true 
height  analysis  to  lower  altitudes  or  lower  electron  densities 
than  possible  from  the  conventional  analyses  alone. 

The  data  were  obtained  utilizing  ionosondes,  a  tnree 
frequency  HF  propagation  measurement  of  both  phase  and  amplitude 
and  riometer  measurements  of  cosmic  noise  absorption.  To  develop 
the  profiles  presented  herein,  the  group  height  records  from 
ionograms  and  phase  height  measurements  from  the  HF  experiment 
were  used. 
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ANALYSIS  OF  LOW  LATITUDE  P  AND  E  REGION  DATA 


I .  INTRODUCTION 


The  analysis  of  high  frequency  phase  and  amplitude  measure¬ 
ments  in  conjunction  with  ionogram  group  path  and  true  height 
electron  density  profile  analyses  was  continued  for  the  data 
obtained  at  Tongatapu  and  Tutuila  in  the  spring  and  fall  of  1962. 

The  goal  of  the  analysis  is  to  better  define  the  electron  density 
profile  at  low  altitudes,  at  electron  densities  and  altitudes  at 
which  the  conventional  true  height  electron  density  profile 
analyses  are  inadequate,  and  to  determine  the  short  term  variations 
in  the  density  profile  in  the  80  to  120  km  region  of  the  ionosphere. 
The  latter  data  is  intended  to  furnish  information  on  the  time 
scale  of  short  term  phase  and  amplitude  disturbances  originating 
in  this  region.  In  order  to  define  short  term  disturbances,  the 
average  profiles  must  be  determinable.  This  effort  comprised 
the  first  phase  of  the  research  effort,  which  extended  over  the 
first  year  of  the  contract. 

Analysis  of  high  frequency  phase  and  amplitude  data  and 
ionograms  in  terms  of  phase  path,  integrated  electron  content, 
and  group  paths  which  are  related  to  an  identifiable  propagation 
mode  (e£  1  hop  F  mode,  0  or  X  wave)  was  initiated  in  the  first 
semi-annual  reporting  period.  The  methods  of  analysis  are 
outlined  in  the  first  semi-annual  report  on  the  contract.  The 
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research  effort  in  the  second  semi-annual  reporting  period 
consisted  of  the  selection  and  analysis  of  data  from  a  typical 
five  day  period,  the  dawn  (1600  to  2100  GMT,  0500  to  1000  LT) 
for  the  days  Octob  r  15  to  19,  1962.  Because  of  the  bulk  of  the 
analyzed  data,  only  a  small  selection  is  presented  herein.  The 
true  height  curves  for  this  data  will  be  presented  in  the  next 
semi-annual  report.  A  shorter  period  of  time,  1700  to  1800  GMT 
on  October  18th,  was  selected  to  represent  the  progress  in  true 
height  analysis  for  a  case  of  one  of  two  solar  flares  which  was 
observable  from  the  low  latitude  sites.  The  profile  changes  and 
gross  effects  on  the  individual  data  were  clearly  discernible 
from  background;  hence  an  interesting  test  (and  somewhat  applicable 
sample  perturbation)  is  obtained  for  the  analysis. 

A  number  of  electron  density  profiles  are  presented  for 
the  period  immediately  after  the  flare.  A  variation  of  some  of 
the  parameters  of  the  analysis  demonstrates,  but  does  not 
quantitatively  evaluate,  the  range  of  accuracy  of  the  analysis. 

In  order  to  evaluate  quantitatively  the  analysis  of  electron 
density  profiles,  it  will  be  necessary  to  compare  electron  density 
profiles  computed  from  ionospheric  models,  and  those  computed  in 
this  analysis  with  in-situ  measurements.  Work  will  continue  in 
the  next  semi-annual  contract  period  to  determine  more  electron 
density  profiles,  and  to  compare  them  with  independently  obtained 
results  such  as  the  ionogram  true  height  analyses  and  other 
input . 
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Upon  initiation  of  the  second  phase  effort,  the  description 
of  short  term  variations  in  the  electron  content  and  profile  in 
the  lower  ionosphere  has  taken  place,  in  that  short  term 
(minute  scale)  variations  which  are  sr.own  in  data  have  been 
reduced  More  data  analyses  in  the  next  reporting  period  will 
show  the  statistical  characteristics  of  the  phase  and  amplitude 
changes  and  hence  the  electron  content  or  profile  variations. 
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II. 


DATA  ANALYSIS 


A. _ F-Region  True  Height  Analysis 

In  order  to  provide  detailed  confirmation  of  an  ionospheric 
model  a  time  series  of  ionograms  was  analyzed  for  the  dawn  period 
(1700-1800  GMT)  on  10  October  1962  using  the  method  of  Paul  and 
Wright. ^  These  analyses  were  carried  out  at  ITSA  by  the  Wright 
group.  Figure  1  illustrates  the  series  of  true  height  profiles 
obtained  for  this  period,  and  Figure  2  illustrates  some  of 
those  obtained  for  the  dawn  period  on  16  October  1962  at  the 
same  site,  Tutuila. 

The  electron  density  profile  above  the  F-region  peak  was 

(2) 

extended  as  described  in  the  first  semi-annual  technical  report'  ' 
by  adoption  of  a  constant  scale  height  Chapman  layer  fit  to  the 
peak  of  the  true  height  analysis  curves. 

The  absorption  variation  for  the  dawn  period  of  18  October 
1962  was  computed  using  the  formulation  d ascribed  in  the  previous 
report,  for  which 

Abs  (30  MHz,  db/km)  =  4.6  x  104  N0  j2  +  V^2)  ; 

the  temperature  dependent  electron-ion  collision  frequency  is 
given  by 


^ei  =  1.8  Ne  Te"3/2  In  (1.54  x  108  Te3/Ne) 
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FIG.  I  TRUE  HEIGHT  ANALYSES  FOR  18  OCTOBER  1962 


30C 


In  this  case  the  electron  temperature  is  assumed  to  be  equal  to 
the  ambient  temperature,  in  which  case  the  computed  absorption 
is  a  maximum. 

B. _ Determination  of  Group  Path  from  Ionograms 

Ionograms  from  the  dawn  period,  1600  to  2100  GMT  lor 
the  interval  15  to  19  October  1962  were  reduced  to  a  virtual 
oblique  group  path  format  by  the  technique  of  Smith^  described 
in  the  first  semi-annual  report.  This  reduction  was  carried  out 
for  the  F  propagated  modes  of  both  magneto-ionic  ordinary  and 
extraordinary  and  high  and  low  ray  cases.  The  frequencies  used 
were  4,  6  and  9  MHz  and  the  propagation  path  was  from  Tonga tapu 
to  Tutuila.  In  order  to  best  describe  the  growth  of  the  E-region 
and  upper  D-region  during  the  dawn  pexiod,  only  the  F  modes 
were  selected  since  they  pass  through  the  lower  ionosphere  largely 
in  a  non-deviative  manner. 

The  results  of  the  five-day  dawn  effect  reduction  effort 
are  illustrated  by  Figure  3,  which  comprises  the  data  for 
19  October  1962.  The  oblique  grcup  path  is  presented  in  kilometers 
versus  time  tor  each  mode  and  for  each  ray  path  when  appropriate. 
Two  distinct  modes  of  behavior  are  present.  The  ray  path  often 
is  not  established  on  either  ordinary  or  extraordinary  traces 
on  the  higher  frequencies.  Therefore,  as  soon  as  the  path  is 
established,  the  virtual  path  length  decreases  rapidly  as  the 
reflection  point  is  lowered  on  the  developing  F-region  density 
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1800 

GMT  TIME,  HOURS 


O-WAVE  GROUP  PATHS,  F-MODE 
1600  -  2100  GMT 
19  OCTOBER  1962 


profile.  Sometime  after  the  dawn  at  E-region  altitudes,  the 
virtual  height  on  all  paths  is  raised  once  more  as  the  increasing 
ionization  causes  retardation  of  the  waves.  This  effect  is 
inversely  proportional  to  frequency  and  the  lower  frequency 
F-modes  are  lost  first. 

In  order  to  confirm  the  ionospheric  electron  density 
profile  model  by  comparison  with  detailed  true  height  curves 
obtained  from  ionogram  analysis,  a  short  period  at  dawn  on  the 
18th  of  October  was  chosen  for  a  detailed  intensive  time  study. 
The  group  paths  for  this  period  are  shown  in  Figure  4. 

In  the  case  of  both  the  longer  term  dawn  effects  and  the 
intensive  study,  the  virtual  group  paths  are  converted  into 
apparent  Doppler  frequency  formats  so  that  they  may  be  directly 
compared  with  the  HF  data  in  the  analysis.  In  addition,  the 
Doppler  is  integrated  from  a  given  starting  time  in  order  to 
give  the  relative  phase  path  change  from  that  time;  this  process 
also  tends  to  smooth  the  otherwise  large  short  term  variations 
in  the  Doppler  caused  by  the  +  2  km  error  of  reduction  of  the 
ionogram  records. 

The  results  of  this  conversion  are  illustrated  in 
Figure  5  for  the  19  October  dawn  study  data  and  in  Figure  6 
for  the  detailed  analysis  of  the  data  for  the  18th. 

C_. _ High  Frequency  Phase  Paths 

The  description  of  the  separation  of  the  modes  of 

propagation  into  F-region  and  E-region  modes  as  outlined  in  the 

1 1 T  RESEARCH  INSTITUTE 


9 


GMT  TIME,  HOURS 


FIG.  4  OBLIQUE  GROUP  PATHS,  F-MO0E 
TUTUILA  IONOGR  AMS 
1700-1800  GMT 
18  OCTOBER  1962 


EGRATED  GROUP  PATHS,  F-MODE 


RELATIVE  INTEGRATED  GROUP  Pi 


previous  semi-annual  report  was  followed  in  order  to  determine 
the  changes  in  the  F-modo  phase  path  lengths  for  the  dawn  periods 
on  15  to  20  October  1962  and,  in  addition,  for  the  detailed 
analysis  period  of  1700-1800  GMT  on  18  October.  Examination  of 
the  phase  records  of  the  HF  data  for  these  periods  showed  no 
unusual  features  in  the  sense  that  the  F-region  wave  could  not 
be  identified.  Separation  of  the  ordinary  and  extraordinary 
waves  could  be  accomplished  by  reference  to  the  oblique  group 
height  records.  That  is,  the  F,  1-hop  mode  was,  in  general, 
established  for  the  extraordinary  wave  before  it  was  established 
for  the  ordinary  wave.  One  can  easily  see  this  behavior  by 
referring  to  Figure  4,  which  illustrates  the  establishment  of 
the  9MHz  extraordinary  wave  group  path  minutes  before  the  ordinary 
wave  group  path  is  established.  The  ionogram  records  and  the 
group  path  data  which  is  obtained  thus  provide  guidance  in 
interpretation  of  the  phase  path  data  in  terms  of  the  particular 
magneto-ionic  mode  being  received.  Mode  changes  and  the  amplitude 
behavior  in  terms  of  p’  3  '  establishment  will  be  discussed  in  the 
next  section. 

Phase  path  data  on  the  three  HF  frequencies,  4,  6  and  9  MHz 
were  reduced  for  the  dawn  period  1600-2100  GMT  and  in  addition, 
a  detailed  minute  by  minute  reduction,  was  carried  out  for  the 
1700-1800  interval  on  the  18th.  These  data  are  illustrated  bY 
Figure  7.  For  example,  identification  of  the  modes 
of  propagation  and  the  continuity  of  the  modes, when  switching 
between  E  and  F-modes  is  taking  place,  is  most  easily  done  when 
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DOPPLER  FREQUENCY,  CPS 


GMT  TIME  ,  HOURS 

FIG.  7  DOPPLER  FREQUENCY,  E-MODE 

1700-  1800  GMT 
18  OCTOBER  1962 
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the  data  is  presented  as  a  Doppler  frequency  instead  of  an 
absolute  phase.  When  the  modes  are  interrupted,  the  Doppler 
frequency  may  be  interpolated  more  easily  than  the  absolute  phase 
however,  the  absolute  phase  record  will  have  an  ambiguity  of 
2  N  71  radians  if  the  record  is  interrupted.  However  the  value 
of  N  is  usually  small  if  the  interruption  is  not  long.  After 
mode  identification  is  established  and  the  necessary  Doppler 
interpolations  are  made,  the  data  may  be  reintegrated  and  the 
relative  changes  in  the  phase  path  of  the  interval  of  interest 
may  be  determined.  This  is  shown  in  Figure  8  with  respect  to 
the  data  for  the  dawn  (1700-1800)  of  the  18th.  At  this  point 
it  is  worthwhile  to  note  that  no  significant  periods  of  loss 
of  a  mode  were  present  in  this  data;  therefore  no  appreciable 
error  is  introduced  through  interpolation  of  Doppler  frequency 
records  through  missing  data  periods. 

D. _ HF  Amplitude  Data 

The  amplitudes  recorded  on  the  high  frequency  receivers 
are  typically  the  sum  of  a  number  of  mode  amplitudes,  the 
ordinary  and  extraordinary  waves  of  the  1  hop  F-mode  plus  the 
ordinary  wave  of  the  E-mode.  The  combination  of  modes  is 
generally  not  so  great  in  number  that  a  random  type  of  fading 
pattern  is  produced.  However,  careful  separation  of  the  effect 
of  the  individual  modes  is  necessary  in  order  to  determine  the 

I  a  T  RESEARCH  INSTITUTE 


15 


amplitude  of  a  single  mode,  which  is  the  quantity  required  for 
electron  density  profile.  The  fading  patterns  typically  consist 
of  two  components:  the  intermode  fading  which  has  a  fading  rate 
close  to  the  Doppler  frequency  of  che  more  rapidly  changing 
mode,  and  the  Faraday  rotation  fading  which  usually  has  a 
frequency  much  lower  than  the  Doppler  frequency  cf  any  of  the 
modes  present. 

In  general,  it  has  been  determined  that  the  presence  of 
strong  E  modes,  wherein  the  Doppler  record  spends  most  of  the 
time  tracing  the  E-modi  signal,  is  exclusive  to  the  case  wherein 
Faraday  rotation  fading  occurs  In  the  F-mode.  In  a  practical 
sense,  this  is  an  important  ex.  lusion  since  the  fading  rates 
for  E-mode  interference  may  be  in  ''he  same  frequency  range  as  those 
expected  from  Faraday  rotation  in  the  F-mode  signal.  For  most 
cases  mentioned,  it  will  be  assumed  thau  che  E-mode  extra¬ 
ordinary  wave  is  of  much  lower  amplitude  than  the  ordinary  wave, 
and  thus  will  not  contribute  appreciably  to  tne  fading  pattern 
or  to  the  overall  signal  strength. 

Another  assumption  which  may  be  shown  qualitatively  is 

that  the  Faraday  fading  depth  will  be  controlled  by  the  presence 

of  the  extraneous  modes  rather  than  by  the  antenna  pattern. 

That  is,  a  dipole  antenna  pattern,  even  over  poor  ground,  has  at 

least  a  10  db  null  for  cross  polarized  waves.  The  Faraday  fading 

signal,  on  the  other  hand,  will  not  be  a  strictly  linearly 

polarized  wave,  although  it  is  the  ^ajor  axis  of  the  ellipse 

of  polarization.  Moreover,  an  important  second  mode  will  usually 
be  present. 
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With  the  aforementioned  assumptions  in  mind,  the  analysis 


of  multi-mode  fading  when  both  Faraday  rotation  and  intermode 
fading  is  present  may  be  illustrated  as  in  Figure  9.  The 
four  parameters  for  this  particular  fading  pattern  through 
are  related  to  the  three  mode  amplitudes  by  the  equations: 

A1  =  Sfo  +  Sf x  +  Se 


A2  =  Sfo  "  Sf x  +  Se 


A3  "  Sfo  +  Sfx  “  Se 


A4  =  Sfo  '  Sf x  "  Se 


In  many  cases  the  intermode  fading  pattern  will  be  deeper  than 
the  Faraday  fading  pattern;  that  is,  the  major  amplitude  difference 
will  be  between  A^  and  A^  rather  than  in  the  envelope  described 
by  A-^  and  A2.  However,  the  method  of  analysis  will  remain  the 
same . 

The  results  of  the  combined  Faraday  and  intermode  fading 
analysis  are  shown  in  Figure  10,  which  was  given  in  the  first 
semi-annual  report.  The  values  of  mode  strength  will  be  used 
later  in  the  analysis  of  the  electron  density  profile  for  the 
period  1700- ^00  GMT  on  October  18th. 
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FIG.  9  FADING  PATTERN  FOR  SIMULTANEOUS  FARADAY  FADING 
AND  INTERMODE  INTERFERENCE  FADING. 
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Total  Signal  -  Minute  Averages 


E-Region  Mode 


III.  THEORY 


A_. _ Theory  -  Group  and  Phase  Paths  in  the  Parabolic  Layer 

A  close  approximation  to  the  fully  developed  Chapman 
layer  (x,  z  0)  is  the  parabolic  layer  described  by  the  equation 

N  =  Nq  (1  -  z2) 

where  z  is  the  reduced  height  =(h  -  h^Y;  h  is  the  actual  height, 
hQ  is  the  height  of  maximum  density,  Nq  is  the  maximum  layer 
density,  and  Y  is  the  layer  semithickness  =  2H.  Rewritten  in 
terms  of  frequency,  the  expression  is  fR  =  fp  (l-(z/2)  ),  where 
fn  is  the  local  plasma  frequency  and  fp  is  the  penetration 
frequency.  The  group  and  phase  paths  for  reflection  of  a  high 
frequency  signal  within  such  a  layer  have  been  shown  to  be 
(Davies,  1965): 

G  -  hQ  -  Y  +  1/2  Y(f/£p)  In  F(fp  +  f)/(fp  -  f) 

and 

P  -  hD  -  (1/2)  Y  -  d/4)  Y  (£p/£  +  £/£p)  In  |(fp  +  9/£p  -  f) 

•  m 

Of  particular  interest  to  this  analysis  is  the  separation 
of  group  and  phase  path  effects  which  occur  in  the  deviative 

1 1 T  RESEARCH  INSTITUTE 


21 


region,  the  parabolic  layer,  from  the  group  and  phase  path 
effects  which  occur  in  the  lower  ionosphere.  The  effects  may 
be  related  to  the  basic  format  of  data  which  is  used  elsewhere 
in  the  report,  namely  the  difference  between  group  and  phase 
paths,  G-P.  This  quantity  is  found  from  the  above  equations 
and  is  G-P  *  Y/4.  Figure  11  illustrates  the  variation  of 
the  reduced  path.,  (G-P) /(Y/4),  as  a  function  of  th>  g»  ment 
f/fp.  It  is  clear  from  the  figure  that  time  changes  in  the 
quantity  G-P,  that  is,  the  measured  value  of  G-P  at  different 
times  will  be  sensitive  to  both  Y  and  f/f^  changes.  These  data 
will  be  obtained  from  ionogram  analysis. 

A  further  element  of  the  use  of  this  theory  in  the 
analysis  is  the  implicit  assumption  of  the  validity  of  the 
equivalent  group  path  theory  of  Martyn  concerning  the  phase 
path.  It  may  be  argued  that  the  phase  path  changes  in  the 
deviative  region  are  a  small  part  of  the  total  phase  path  change. 
Thus  use  of  the  equivalent  vertical  frequency  as  in  the  group 
path  analysis  will  be  allowable  The  error  induced  by  this 
assumption  will  be  investigated  further  and  an  attempt  will  be 
made  to  prove  an  equivalent  statement  for  phase  path. 

_ Group  and  Phase  Paths  in  the  Exponential  Layer 

The  lower,  non-deviating  region  of  the  ionosphere  will 
be  approximated  by  an  exponential  function, 

fj  =  F2  exp  (az)  , 
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FIG. 11  VARIATION  OF  GROUP-PHASE  HEIGHTS  VS 

FREQUENCY  FOR  PARABOLIC  LAYER  REFLECTION. 


1 
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where  z  is  again  a  reduced  height,  F  is  the  plasma  frequency  at 
a  reference  height,  and  f  is  the  plasma  frequency  at  any  height 
in  the  layer.  In  all  cases  in  the  analysis  the  transmitted 
frequency  f  will  be  much  larger  than  the  highest  plasma  frequency 
in  the  layer  so  that  the  layer  may  be  treated  as  a  strictly 
non-deviating  region.  Therefore,  for  this  case  the  group  and 
phase  paths  are 

G  -  2/a  In  j^f/F  +  {  (f/F)2  -  1  }  - 
P  -  G  -  (2/a)  (F/ f )  {(f/F)2  -  l}  * 
which  become  in  the  limiting  case  f>}  F, 

G  =  L 

P  -  L  ’  2/a  . 

As  before,  the  difference  in  phase  and  group  paths  may  be 
determined  as 


G-P  =  (2/a) (F/f )  {(f/F)2  -  l}  *  . 

This  quantity  is  plotted  in  Figure  12  as  a  function  of  reduced 
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GROUP  AND  PHASE  PATH  VS  FREQUENCY 
flAL  LAYER  TRANSMISSION. 


variables  (G-P)/(2/a)  and  f/F,  In  this  case,  f/F  >  1  is 
always  true. 

Of  more  importance,  the  actual  reduced  phase  path  for 
the  layer  is  also  shown  in  Figure  12.  It  may  easily  be  seen 
that  for  an  equivalent  phase  path  change,  the  group  path  change 
will  be  a  constant  factor  over  much  of  the  frequency  range  of 
interest.  The  magnitude  of  the  group  path  change  will  be  larger 
than  the  phase  path  change  for  f/F  3.6;  because  of  the  lack 
of  sensitivity  of  the  i.onogram  records  to  wavelength  scale  change 
of  path,  the  phase  path  will  contribute  all  to  the  measurable 
data.  The  group  and  phase  paths  will  be  assumed  to  be  without 
deviation  over  the  entire  exponential  region  and  the  oblique 
propagation  factor  will  be  taken  into  account  by  suitable 
correction  of  the  height  terms  in  the  formulation, 
ie  h (oblique)  *  hv  sec  j 6. 

C_. _ Explicit  Variations  of  Layer  Parameters 

Insofar  as  the  analysis  is  concerned  with  defining  a 
change  in  E-region  electron  profile  description,  the  parameters 
of  the  description  are  the  time  variables;  that  is,  the  value 
of  a  and  F  will  change  with  time.  Thus  the  sensitivity  of 
the  group  and  phase  paths  to  variations  of  these  parameters 
must  be  determined.  The  following  equations  describe  these 
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variational  terms: 


Parabolic  Layers 
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Exponential  Lavers 
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IV. 


DETERMINATION  OF  ELECTRON  DENSITY  PROFILES 


Phase  path  data  on  4  and  6  MHz  were  used  in  conjunction 
with  the  parametric  description  of  the  parabolic  layer  outlined 
in  the  previous  section  to  determine  the  electron  density  profile 
for  certain  times  between  1720  and  1800  GMT  on  18  October  1962. 
The  data  utilised  for  this  analysis  was  described  in  Section  II 
and  the  true  height  analyses  of  ionograms  taken  at  simultaneous 
times  were  utilized  to  determine  the  plasma  frequency  or  altitude 
of  the  upper  boundary  point  of  the  layer.  In  many  cases,  when 
only  the  observed  E-region  plasma  frequency  was  ured,  the  peak 
height  parameter  hQ  and  the  thickness  parameter  Y  became  too 
large.  This  is  because  the  electron  density  near  the  peak  of 
the  E-region  is  better  represented  by  the  parabolic  layer 
description  than  the  density  near  the  respective  reflection 
points  for  the  4  and  6  MHz  waves.  However  choice  of  a  layer 
height  parameter  hQ,  from  the  true  height  analyses,  allowed 
computation  of  a  realistic  thickness  parameter,  Y,  from  both  sets 
of  data. 

The  electron  density  contours  computed  from  the  parabolic 
layer  approximation  and  application  of  only  the  plasma  frequency 
at  the  peak  of  the  parabola  are  illustrated  in  Figure  13.  It 
may  be  seen  that  agreement  of  the  computed  profiles  from  the  true 
height  profiles  obtained  from  ionograms  becomes  poor  as  the 
computations  are  repeated  and  carried  backwards  in  time  from 
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Che  1759  prof'le.  This  backward  computation  was  carried  out 
instead  of  a  icrward  time  method  because  the  1759  profile 
appeared  to  be  the  best  known.  The  direction  of  the  time  or 
iteration  interval  is  immaterial  since  the  phase  and  amplitude 
data  can  be  carried  either  way  from  a  reference  point.  The 
profile  comparison  made  at  1759  resulted  from  fitting  the  three 
parameter^  hQ  V  and  f  to  the  ionogram  derived  true  height  curve. 

Thus  it  would  be  expected  that  the  thr  e  curves,  1759,  1750  and 
1740  would  agree  quite  well  with  the  ionogram  results.  This  is 
not  the  case  and  as  mentioned  before,  results  from  the  sensitivity 
of  the  computation  to  the  exact  choice  of  hQ. 

Before  moving  on  to  mere  accurate  techniques  of  fitting 
the  profiles,  it  must  be  emphasized  that  the  failure  of  the 
initial  attempt  is  not  in  the  method  but  in  the  extension  of  the 
profile  from  where  it  fits  in  the  4  and  6  MHz  region  to  where 
it  is  forced  to  fit  the  ionogram  profile  results  a*:  higher  altitudes. 
In  each  case  this  region  of  joining  is  appreciably  above  the 
altitude  of  the  4  and  6  MIz  reflection  altitude.  Therefor’, 
if  the  contour  as  initially  derived  is  considered  to  be  that 
applicable  mainly  to  the  4  to  6  MHz  region,  it  will  be  quite 
representative.  The  goal,  however,  is  to  produce  a  continuous 
contour  which  will  extend  the  ionogram  derived  density  contour 
u..wrwards  in  an  accurate  fashion. 

In  order  to  better  fit  the  results  of  the  hf  analysis 

with  the  true  height  contour,  the  parameter  hy  is  varied  and  the 

profile  fit  is  examined.  Figure  14  shows  the  effect  of  vaiiation 
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of  hQ  on  the  1720  GMT  profile.  It  may  be  seen  that  the  layer 
thickness  which  is  most  comparable  to  that  expected  from  the 
other  data  also  corresponds  to  the  best  fit  of  the  true  height 
curve . 

In  conclusion,  since  the  desired  results  are  the  description 
of  the  ionosphere  at  heights  or  densities  below  those  amenable 
to  true  height  analysis,  the  best  ionogram  and  hf  true  height 
curves  are  presented  in  Figure  15. 

One  qualifying  note  must  be  added;  the  matching  region  of 
the  hf  and  ionogram  true  height  curves  depend  upon  the  accuracy 
of  the  ionogram  true  height  analyses  for  the  initial  or  starting 
point  of  the  analysis.  Thus,  the  quality  01  the  hf  true  height 
profiles  will  depend  directly  upon  the  quality  of  the  ionogram 
true  height  profiles  and  hence  will  possess  some  degree  of 
ambiguity.  However,  use  of  the  amplitude  data  will  impose 
additional  altitude  constraints  on  the  profiles  and  will  aid  in 
fixing  lower  boundary  conditions  just  as  the  ionogram  results 
are  used  for  upper  boundary  conditions.  It  was  impossible  to 
apply  this  extended  analysis  to  this  time  period  because  no 
4  or  6  MHz  amplitude  data  were  obtained  because  of  equipment 
difficulty. 
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IONOGR AM  TRUE  HEIGHT  CURVE 
1725  GMT  18  OCTOBER  1962 


V. 


IONOSPHERIC  DISTURBANCE  CHARACTERIZATION  FOR  ALL  DATA 


PERIODS 


Data  records  from  two  intervals  were  obtained  in  the 
South  Pacific  sites  during  the  1962  test  series.  The  first 
period  was  in  the  spring  and  good  data  was  obtained  from 
July  4  to  11  and  from  July  21  to  26.  In  the  fall,  a  few  days 
data  were  obtained  in  September;  however  most  of  the  useable 
data  covered  the  period  15  October  to  5  November.  The  data 
chosen  for  initial  examination  for  this  analysis  was  in  the 
15  October  to  20  October  interval. 

It  is  desirable  to  characterize  the  degree  of  solar, 
geomagnetic  and  ionospheric  disturbance  for  the  days  in  which 
data  were  obtained  so  that  a  representative  choice  of  analysis 
periods  can  be  made.  In  particular,  the  differences  in  the 
statistical,  or  short  term,  phase  and  amplitude  variations  in  the 
HF  data  between  high  ionospheric  disturbance  days  and  normal  or 
quiet  days  is  an  important  factor  to  bi_  examined  during  the 
second  phase  of  the  study. 

In  order  to  allow  this  selection,  the  worldwide  geomagnetic 

disturbance  index  K  ,  the  solar  radio  noise  index  S,  and  the 

P 

existence  of  solar  flare  activity  are  tabulated  in  this  section. 

In  addition,  the  acquisition  index  numbers  of  the  appropriate 
data  are  listed  in  order  to  allow  ease  in  locating  particular 
samples  of  data. 
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Examination  of  the  tables  will  indicate  that  the  period 
of  detailed  ionospheric  profile  analysis  (1700-1800  GMT,  101862) 
corresponds  to  one  of  the  two  opportunities  present  during  the 
entire  experiment  to  measure  the  ionization  produced  by  a  solar 
flare.  Unfortunately  the  flare  took  place  during  early  morning 
and  thus  the  full  effects  at  lower  altitudes  were  probably  not 
realized  because  of  the  grazing  incidence  of  the  x  radiation. 
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TABLE  A 


SOLAR  DISTURBANCE  INDEX  GEOMAGNETIC  DISTURBANCE  INDf.X 
FOR  JULY  AND  DATA  REFERENCE 


Date 

S(1) 

Sum  K(J( 

2)  Tut . ABS 

Tonga  ABS 

11 F6 

liFb 

11 F  9 

1 

18 

1118 

15q 

1118 

3 

r 

5 

90 

88 

16q 

251) 

26  D 

1118 

1.118 

1118 

1189 

1188,93 

1187,91 

1192 

1186 

6 

86 

21 

11 18 

1185 

7 

88 

15 

1118 

1195 

1196 

1 190 

8 

83 

20 

1118 

1198 

1200 

1202 

9 

80 

12Q 

1118 

1197 

1196 

1203 

10 

81 

L  J 

1118 

1199 

1201 

1206 

11 

83 

18 

1160-1 

1161 

1206 

1205 

1207 

12 

82 

15q 

1651 

1161 

13 

86 

19 

1651 

1161 

16 

86 

17 

1651 

1161 

15 

85 

llq 

1651 

1161 

lb 

86 

6Q 

1651 

1161 

17 

86 

6Q 

1651 

1161 

18 

82 

8(7 

1651 

1161 

19 

80 

18 

1651 

20 

80 

26 

1651 

21 

79 

21 

1651 

1519 

1518 

1517 

22 

80 

16q 

1651 

1519 

1518 

1517 

23 

78 

16 

1651 

1612 

1523 

15  22 

1521 

26 

78 

20 

1651 

1613 

1527 

i '2*: 

1525 

25 

76 

20 

1651 

1613 

1527 

1526 

1525 

26 

76 

37D 

1122 

1613 

1531 

1530 

1529 

27 

76 

3  ID 

1617 

28 

76 

26D 

1617 

29 

73 

17 

1625 

3<i 

72 

10Q 

3  i 

7  3 

15 

Data  is  daily  values 
taken  from  "Handbook 
McCraw  Hill,  1965. 

of  solar  flux  at 
of  Geophysics  anc 

2800  MHz  in  units  of 
Space  Environments" 

.n 

10  ““Win 
ed.  S . 

2hz_1 

Valley 

2.  Data  is  Sum  K..  index  value  Cron  J.  V.  Lincoln,  j.  Geophys.  Res.  67,  6875, 
1962.  Q  and  D  indicate  5  selected  quiet  or  disturbed  days,  q  antT"d 
indicate  10. 

J.  Accession  numbers  correspond  to  chose  in  the  ionospheric  data  catalog 
at  IIT  Research  Institute. 
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TABLE  R 


SOLAR  AND  GEOMAGNETIC  DISTURBANCE  INDICES  FOR 
OCTOBER-NOVFMBER  DATA 


Date 

s<i) 

Sum  Kp(2) 

Tut.  ABS 

Tonga  ABS 

HF4 

HF6 

HF9 

October 

13 

95 

18q 

2228 

1576 

2269 

2270 

2268 

14 

95 

33 

2228 

1576 

2269 

2270 

2268 

15 

94 

18Q 

2228 

1576 

2272 

2273 

2274 

16 

91 

25 

2228 

1588 

2288 

2287 

2286 

17 

91 

13Q 

2228 

1588 

2275,92 

2276,91 

2277,90 

18 

89 

22 

2387 

1588 

2280,96 

2278,95 

2279,96 

19 

87 

27 

2387 

1588 

2280 

2278 

2279 

20 

88 

19q 

2387 

1588 

2284 

2283 

2282 

21 

87 

22 

2387 

1588 

1914 

1913 

1912 

22 

85 

29 

1656 

1588 

1914 

1913 

1912 

23 

84 

29 

1656 

1588 

1919 

1916 

1917 

24 

89 

31 

1656 

1977 

1922 

1921 

1920 

25 

87 

35D 

1656 

19^ 

1 Q22 

1921 

1920 

26 

87 

34D 

1762 

J.9/  / 

1842 

1850 

1841 

27 

86 

31 

1762 

1884 

1845 

1844 

1843 

28 

82 

26 

1762 

1884 

1845 

1844 

1843 

29 

80 

22 

1762 

1884 

1849 

1848 

1847 

30 

82 

23 

1762 

1884 

1354 

1853 

1852 

31 

81 

loQ 

1762 

1884 

1858 

1857 

1856 

November 

1 

80 

16 

1762 

1888 

1866 

1865 

1864 

2 

80 

22 

1762 

1894 

1893 

1892 

3 

90 

23 

.1659 

2157 

1894 

1893 

1892 

4 

82 

26 

1659 

2157 

1901 

1900 

1899 

5 

82 

llq 

1659 

2157 

1901 

1900 

1899 

1.  Data  Is  daily  values  of  solar  flux  at  2800  MHz  in  units  of  10  ^  Wm  ^hz  ^ 
taken  from  "Handbook  of  Geophysics  and  Space  Environments"  ed.  S.L. Valley, 
McGraw  Hill,  1965. 

2.  Data  is  Sum  Kp  index  value  from  J.V.  Lincoln,  J.  Geophys.  Res.  67,  4875, 
1962.  Q  and  D  indicate  5  selected  quiet  or  disturbed  days,  q  ancT  d  indicate 
10. 

3.  Accession  numbers  correspond  to  those  in  the  ionospheric  data  catalog 
at  IIT  Research  Institute. 
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TABLE  C 

SUMMARY  OF  SOLAR  FLARE  ACTIVITY 


Date 

Beginning  Time 

Ending  Time 

Importance 

July 

5* 

1937Z 

1947Z 

2 

8 

1540Z 

1604Z 

1+ 

October 

15 

1525Z 

1550Z 

1+ 

16 

1000Z 

1038Z 

1+ 

* 

17 

0717Z 

Not  Obs. 

1+ 

* 

18 

1705Z 

1730Z 

20 

1625Z 

Not  Obs. 

1+ 

21 

1535Z 

1645Z 

1+ 

23 

1705Z 

Not  Obs. 

1+ 

23 

1642Z 

1745Z 

2 

24 

1306Z 

1346Z 

End  Not  Obs. 

1+ 

28 

1751Z 

Microflare 

Indicates  flares  which  could  be  observed  from  the  Pacific  low  latitude  sites 


UPPER  PROFILE  FROM  IONOGR AM  ANALYSIS 
LOWER  PROFILE  FROM  HF  PHASE  PATH  ANALYSIS 


FIG.  15  ELECTRON  DENSITY  PROFILES  FOR  18  OCTOBER  1962,  1725  TO  1800  GMT. 
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APPENDIX  A 


ERROR  ANALYSIS  OF  GROUP  AND  PHASE  PATH 
ANALYTICAL  APPROXIMATIONS 


In  analysis  of  the  down  increase  in  ionospheric  electron 
content,  the  quantity 

p-g-  /;  l.  -m)  d  P 

j  ij 

is  derived  from  experimental  data.  The  gj.oup  path  is  obtained 
by  reduction  of  ionogram  records  and  the  phase  path  is  determined 
by  hf  phase  measurements.  In  order  to  determine  the  possible 
error  in  excess  phase  path  change,  i.e.  that  part  of  the  phase 
change  determined  by  the  lower  ionospheric  electron  density 
profile  changes,  the  error  in  group  path  measurements,  phase 
path  measurements  and  the  analytical  theory  must  be  determined. 

Measurement  errors  and  the  error  induced  by  data  reduction 
techniques  in  the  virtual  group  path  data,G  ,  are  estimated  to 
be&G  ■  +  2  km.  This  is  equivalent  to  +  26\  ,  +  40X  ,  and 
+  60\  ,  at  4,  6  and  9  MHz  respectively.  This  error  includes 
both  the  geometrical  path  error  and  that  induced  by  wave 
retardation.  The  down  period  F  mode  is  characterizea  by  many 
cycles  phase  change  so  with  a  phase  error  of&P  +1  cycle, 

&  P/p  <  1/p,  and  thus  the  relative  error  is  small  compared  with 
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the  LG/G  +  2/r,  value  as  Tong  as  P  is  great  enough.  For  this 
reason,  the  time  intervals  covered  by  the  measurements  are  made 
large  enough  to  ensure  that  J./P  and  2/G  are  small. 

Errors  may  be  induced  by  the  approximate  expression  chosen 

1/2 

for  the  value  of  w  =  (1  -  X)  '  .  For  example,  if  the  non  deviative 
approximation  is  adopted,  u  =  (1  =  1/2  X),  and  the  expression  mav 
be  in  error  in  the  reflection  region. 

The  quasi- longitudinal  case  is  studied  herein  by  expanding 
1/2 

the  function  ^  =  (1-X)  '  to  the  third  power,  then  evaluating  a 
perturbation  imposed  at  the  largest  usual  value  of  X.  That  is, 
the  expression  (u  +  £u)  is  evaluated  at  X  =  1/2,  its  larges." 
nominal  value  for  F-mode  reflection.  Thus  u(min)  ^.0.7.  The 
ionogram  true  height  analyses  were  evaluated  for  all  three 
frequencies  in  the  period  1700-1800  on  October  18  and  others 
which  were  available  were  examined.  In  most  cases  the  value 
X  =  f0/F  <  0.5. 

The  expansion  of  is: 

'  0 

u  -  Ho  +  Uo  +1  n  "  X 

T: 

1  12  11 

U  =  1  -  X  +  -g  X  -  X  etc. 


u  1st  order 
a  true 


0.75 

0.7 


1  r\  -7 


i.  .  U  / 
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Thus,  use  of  the  quasi- longitudinal  non-deviation  index  of 
refracticn  induces  at  most  an  error  of  a  few  percent  in  the 
phase  path  relationship  to  integrated  electron  content. 

The  group  path  error  is  also  of  interest  but  will  not 
enter  strongly  into  the  error  of  the  analysis  because  the  group 
path  effect  is  subtracted  from  the  data  and  is  not  used 
specifically  in  analytical  form 


—  =  (1  -  X)  ~1^2  =  1  +  7  x  +  | 


15 

8  V  3! 


—  at  X  —  —  1.42  and 

u  i 

u(lst  order)  _  1.25  =  x  12 

,  (true)  1,4 

A  12  percent  induced  error  is  present  in  the  first  order 
approximation  to  the  group  path  analysis  if  group  path  data 
were  to  be  analyzed  in  this  fashion. 
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height  measurements  from  the  HF  experiment  wire  used. 
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